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ABSTRACT 

We describe a new method for simulating ionizing radiation and supernova feedback in galaxy simulations. 
In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique 
to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the 
fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self- 
regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from 
star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. 
By simulating a galactic halo of 2.3 x IO^'M©, we find that the galactic escape fraction, /esc, fluctuates between 
0.08% and 5.9% during a ^ 20 Myr period with a mean value of 1.1%. The flux of escaped photons is not 
strongly beamed, but manifests a large opening angle of more than 60° from the galactic pole. Further, we 
investigate the escape fraction per SFMC particle, /esc(0' ™d how it evolves as the particle ages. We discover 
that the galactic escape fraction is dominated by a small number of SFMC particles with high /esc(0- 
average, the escape fraction from a SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle 
lifetime, 6 Myrs. This is because SFMC particles drift away from the dense gas clumps in which they were 
born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova 
feedback. The framework established in this study brings deeper insight into the physics of photon escape 
fraction from an individual star-forming clump, and from a galaxy. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: star clusters — galaxies: starburst 
— stars: formation — stars: evolution — ISM: HII regions 



1. INTRODUCTION 

Lying at the center of a web of processes that drive cosmic 
evolution is galactic star formation. Decades of observations 
have yielded a number of empirical rules about how it oper- 
ates. The accumulated observational knowledge considerably 
exceeds our physical understanding, and as a result simula- 
tions and semi-analytic models of galaxy formation, reioniza- 
tion, and interstellar medium (ISM) dynamics often rely on 
subgrid models of uncertain validity for converting gaseous 
material into stars and then injecting energy, momentum, and 
newly-synthesized elements back into the ISM. 

Theoretical understanding of star formation has been ham- 
pered primarily by the fact that star-forming gas clumps do 
not exist in isolation. Rather, they are a part of a larger galac- 
tic structure. Molecular clouds are strongly associated with 
peaks of the atomic hydrogen distribution in galaxies (e.g. 
iHuphes et al.ll20^1 0). and both direct kinematic measurements 
and indirect statistical evidence suggest that giant molecular 
clouds (GMCs) spend much of their lifetime accretin g mass 
from their surrounding atomic envelopes (e.g. Fuku i et aH 
|2009). Gas accretion alters molecular clouds by providing 
mass to replace that ejected by stellar feedback, and by adding 
energy to off set the decay of turbulence and stall gravitational 
collap se (e.g. lVazquez-Semadeni et all2010l : lGoldbaum et alj 
1201 Ih . Consequently, the internal processes that drive and 
regulate star formation in molecular clouds cannot be under- 
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Stood outside the larger galactic context. 

Needless to say, numerical investigations of star-forming 
gas clumps in a galactic context pose a problem of inher- 
ently multi-scale nature. In local galaxie s, GMCs have typical 
surface densities of E lOOMgpc"^ dBolatto et alj r2008). 
and star form ation occurs inside even denser clumps (E ^ 
lOOOM.^pc'^: iPlume et a l. 1997; Lombardi et al. 2010). In 
contrast, GMCs' atomic hydrogen envelopes are less dense 
(E ^ 10 Mfjpc^^), and the surrounding ISM can be even more 
diffuse (E < 1 MqPC^^). Indeed, the large dynamic range of 
the problem have prohibited theorists from numerically re- 
solving individual star-forming events in a simulation of a 
long-term galactic evolution. 

A variety of methods has been proposed to limit the dy- 
namic range of a galactic ISM simulation. Simply speaking, 
these are attempts to model star formation and feedback at the 
numerical resolution one can afford. A popular choice among 
these is to use a subgrid model to estimate the pressuriza- 
tion of gas on large scales (^ 10 pc) and impose a pr essure 
floor at the calculated level (e.g. ISpringel et ani2005f) . One 
may also impose artifi cial pressurization to stabilize th e ISM 
at smaller scales (e.g. iTevssier et al]|2010l : FBournaud et alj 
|2010). While these pressure floor models successfully halt 
gravitational collapse and produce equilibrium galaxies, the 
treatment of subgrid physics is model-dependent and must be 
tuned to match observations. Another option is to allow the 
gas to collapse unimpeded and then assume that the gas that 
exceeds the highest resolvable den sity forms a b ound struc- 
ture (Lagrangian sink particle; e.g. iLi et alJl2005(). However, 
this approach may overestimate the amount of gas in the dense 
ISM, especially in simulations that do not include proper feed- 
back. A related approach is to allow unimpeded collapse but 
not to put all the collapsing gas into a sink particle, and instead 
rely o n star formation recipes (Shettv & Ostriker 2008^ Tasker 
& Tan 120091; iDobbs et al.l 120111). This approach, however. 
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Fig. 1 . — Two dimensional schematic views of stai-forming molecular cloud (SFMC) particle formation and two different channels of feedback. Left: formation 
of SFMC particles desciibed in ^2.21 The numerical resolution of 3.8 pc is in accord with the Jeans length for a gas clump of ;i = 1000 cm^' at ~100 K, at which 
point a molecular cloud particle of 1000 M ;, spawns. Middle: radiation feedback described in i]2.4l Rays of io nizing photons carrying the monochromatic energy 
of 1 6.0 eV are adaptively traced for particles of age T = to 6 Myr. Right: supernova feedback described in i|2.3l Along with the radiation, the thermal energy 
by supemovae explosion is injected into the gas cell in which a SFMC particle of age T = 4 to 6 Myr resides. 



necessarily involves a violation of the Jeans resolution con- 
dition at the fin est level, with unknown consequences (Tru- 
elove et al. 1997|). Hence, one might arguably claim that there 
has been no numerical simulation that can self-consistently 
couple the properties of a galaxy to the internal properties of 
star-forming gas clumps. At the same time, no simulation 
has formed a population of numerically converged molecular 
clouds that agrees with the GMC populations observed in ex- 
ternal galaxies. Because of these outstanding challenges, it is 
unclear how the properties of molecular clouds are imprinted 
with the properties of their host galaxies. 

The goal of this work is to improve this situation by per- 
forming first-principles simulations of the behavior of star- 
forming gas clumps in a galactic context. In particular, we 
aim to address the physics of star-forming clumps in a galaxy 
formation simulation with ever-improving computational res- 
olution. As the numerical resolution of galaxy simulations 
dramatically improves, spatial resolution of ^ pc is now fre- 
quently employed and tested by a number of authors (e.g. 
Tevssier et all 120101: lBournaud et~an 120101: iHopkins et all 
20111) . Yet, the formation and feedback of star clusters in 
galaxy formation simulations have been described only phe- 
nomenologically at best. Such examples include placing a star 
cluster particle using an observed star formation density rela- 
tion (e.g. psFR Pgas; ISchmidtllT959l: iKennicutti [19981) . or 
stopping the cooling near a star cluster particle right after its 
birth in an attempt to mimi c the supernova energy injection 
(e.g. iGovernato et al.ll2012l and references therein). These 
phenomenological prescriptions may have worked effectively 
and efficiently at low resolution. However, a careful investi- 
gation is greatly needed to depict the realistic physics of star- 
forming gas clumps in a simulation that resolves the molec- 
ular cloud scale of Ax ^ pc with a corresponding timestep 
of Af ^ 3 X lO'* yr at 10^ K. One should note that this fine 
timestep now starts to resolve the lifetime of typical molecular 
clouds and massive stars, placing many high-resolution simu- 
lations in a great need of temporally-resolved stellar feedback 
model (f or recent attempts to tackle this challenge, see Hop- 
kins et al. l201 lUStinson et al.ll2012h . 

In light of these developments, we take advantage of the 
high dynamic range made possible by adaptive mesh refine- 
ment (AMR) technique to simulate the behavior of < 10000 
star-forming molecular cloud (SFMC) particles and their em- 



bedding galaxy in a single self-consistent numerical frame- 
work. We have implemented a sophisticated description of 
stellar feedback by combining the ultraviolet radiation from 
star clusters and the thermal energy injection by supernova 
explosion. The numerical resolution of 3.8 pc helps to resolve 
the star-forming environment in a 2.3 x 10 Mq galactic halo, 
and to describe the capture of stellar radiation energy. Our 
work funda mentally differs from the representation of stellar 
feedback in lHopkins et al.l (1201 1) and their subsequent stud- 
ies. Specifically, we do not artificially pressurize the ISM at 
small scales, nor rely on a subgrid model to handle the stel- 
lar radiation feedback. Instead, we actually solve the transfer 
equation on the fly. This means that our approach does not 
rely on free parameters to describe the interaction of radiation 
and the gas, and also that we can address questions that previ- 
ous studies cannot. In this paper we focus on one of these: by 
solving the transfer equation, we are able to study the escape 
of ionizing radiation from galaxies self-consistently, with the 
ionizing radiation actually included in the simulation rather 
than added in post-process. Future work will focus on other 
aspects of our simulations, including regulation of star forma- 
tion by ionizing radiation feedback, and mock observations of 
ionized gas tracers such as Ha emission. 

In a series of papers, we will discuss the state-of-the-art 
model that describes feedback from SFMC particles and how 
it enhances our understanding of galactic star formation by 
significantly expanding the usage of numerical simulations. 
The first paper of the series is organized as follows. The 
physics in the simulation code is described in S|2] and the 
initial condition of the numerical experiments is explained in 
21 Sloverviews the performed simulations of dwarf galaxies 
with radiation feedback, and f|5]is devoted to the results of one 
of our experiments, with a strong emphasis on the escape of 
ionizing photons. Particular attention is paid to the evolution 
of the escape fraction per SFMC particle. Assembled in S|6] 
are the summary and conclusions. 

2. PHYSICS IN THE CODE 

We utilize the Eulerian adaptive mesh refinement Enzo, a 
parallel hybrid N-bodv radiation hvdrody namics code (Bryan 
& Norman 1997: Nor man & Br van 1999: Norman et al. 2007|: 
,Wise & AbeI201 1.) . The code contains several basic modules. 
The hydrodynamics module solves the Euler equations of 
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compressible gas dynamics using either a piecewise-parabolic 
method or a ZEUS scheme. The gravity module computes 
the gravitational potential generated by the fluid, collisionless 
dark matter and stellar particles. The latter two are discretized 
onto the grid using an adaptive particle-mesh method, allow- 
ing Enzo to solve the Poisson's equation for the potential via 
fast Fourier transform and/or multigrid solvers. The radia- 
tion transfer module propagates ionizing photons through a 
spatially-adaptive tree of rays drawn around each source, and 
updates the thermal and chemical states of the gas as a result 
of its interaction with these photons. All these modules oper- 
ate on top of an AMR hierarchy that allows users to specify ar- 
bitrary refinement conditions. In addition to these basic mod- 
ules, our version of Enzo-2.1 ^ contains all releva nt physics 
previously considered in galactic simulations (e.g. lKim et alj 
I2OO9. ,201_1) as well as other enhanced physics discussed in 
detail below. As a result of the wide range of physics in- 
cluded, we are well-poised to investigate the physics of star- 
forming clumps in a galactic context. 

2. 1 . Hydrodynamics, Refinement, and Chemistry 

We use the ZEUS hydrodynamics module included in Enzo 
to solve the Euler equations for the collisional fluid (Stone & 
Norman 1992a. b; Anninos & Norman 1994). Dark matter and 
stars are treated as collisionless particles. The grids are adap- 
tively refined by factors of 2 in each axis on gas and particle 
overdensities. The mass thresholds, Mref, above which a cell 
is refined are functions of a refinement level / as 

j^l ^ ^-0.820/ yy^O .IX 

ref.gas ^ ^"ref.gas ' ' 

j^l _^ -0.533/ yy^O (J. 

^"ref,part ~ ^ ^"ref.part V^/ 

for gas and particles, respectively, with M^^^^^^ — M^^f_p^^ = 

4.87 X 10^ Mq for a 32^ root grid. This makes the simula- 
tion super-Lagrangian, refining grids more on small scales. 
At the finest level / — 13 of 3.8 pc resolution, a cell would 
have been further refined if more than 3000 in gas or 
40000 Mp, ~ 2.5 Mdm in particles. This numerical resolution 
is in accord with the Jeans length for a dense gas clump of 
n — 1000 cm^^ at 100 K, at which point a corresponding 
Jeans mass of 2000 Mq collapses to spawn a SFMC particle 
(see m.2i . Because 1000 to 1500 Mq is instantly removed 
from the cell every time a particle is created, the gas mass in 
the finest cell never reaches the refinement threshold. 

We choose the hydrodynamic Courant-Friedrichs-Lewy 
safety number of 0.3 so that the marching timestep is less 
than the time for the advecting wave to reach the adjacent 
cell. We make use of Enzo's multi-species non-equilibrium 
chemistry module to track the reactions among six species 
(H, H+, He, He+, He++, e^), and also employ its cooling 
module to calculate radiative losses due to collisional exci- 
tation cooling, collisional ionization cooling, recombination 
cooling, Bremsstrahlung cooling, and CMB Co mpton cooling 
for hydrogen and helium (lAnninos et al.|[T997l) . Furthermore, 
cooling by metals is included in gas above 10^ K using the 
cooling rates tabulated by Sutherland & Dopita (1993), and 
in gas below lO'* K using the cooling rate approximation of 
iKovam a & Inut suka (|2002|) (see however Vazquez-Semadeni 
et al. (2010) for an important typo). Photoionization heat- 
ing by the metagalactic ultraviolet background of quasars and 
galaxies is also taken into account without consider ing the ef- 
fect of self-shielding dHaardt & Madaulll996l[200lh . 

' http://enzo-project.org/ 
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Fig. 2. — The e volut ion of SFMC particle mass, Muc, as a function of 
particle age T. See i|2.3l and Eq.(4) for detailed explanations. 

2.2. Formation of Star- forming Molecular Cloud ( SFMC) 
Particles 

We describe the model of formation and feedback of SFMC 
particles in the remainder o f |2l Our SFMC particle formation 
is bas ed on ICen & Ostriker (1992) with modifications simi- 
lar to lKim et al.N2011h . Wifli an efficiency of e* = 0.5, the 
finest cell of size Ax — 3.8 pc and gas density Pgas produces 
a SFMC particle of initial mass M^^q — e^PgasAx^ when (a) 
the proton number density exceeds the threshold nthres — 1000 
cm^^, (b) the velocity flow is converging, (c) the cooling 
time fcooi is shorter than the gas dynamical time t^y„ of the 
cell (^ 1 Myr), and (d) the particle produced has at least 
■'^MC > ■'^thres = IOOOMq. As a result, the gas in the finest 
cell is instantly converted into a particle as soon as the cell 
has accumulated more than Mthres/£* = 2000 M0, the Jeans 
mass at n = 1000 cm^^ at 100 K (see the left panel of Fig- 
ure [T]i. This guarantees that a SFMC particle is created be- 
fore an unresolved gas clump interferes with the consistency 
of hydrodynamics (see m.lj . Because we deposit a particle 
when a gas cell of a typical molecular cloud size actually be- 
comes Jeans unstable, the particle in our simulation represents 
a star-forming molecular cloud that is self-gravitating, is thus 
decoupled from the gas on the grid. Additionally, two SFMC 
particles merge in a way that conserves mass and momentum 
if they are separated by less than 5 pc in space, and 3 x 10^ yr 
in age. This is to reduce the number of active SFMC particles 
and thereby expedite the radiative transfer calculation. 

2.3. Evolution of SFMC Particles and Supernova Feedback 

In our simulation, the evolution of a SFMC particle of mass 
MJJjJJ. describes the following stages: [Stage I] At the birth of 
a SFMC particle, 24% of MJ^ji is instandy turned into stars. 
The rest, 0.76MJJJji, is considered as atomic and molecular 
gas which does not participate in star formation, modeling 
the inefficient star f ormation in molecular c louds observed 
in local clusters (e.g. iKmmholz & Tanll2007l a nd references 
therei n) and predicted by numerical studies (e.g. lMurrav et al.l 
|2010|) . [Stage II] Assuming a Salpeter initial mass function 
(IMF;|SalEete3[i955l) between 0. 1 and 300^^0,*^ 17% of stel- 
lar mass, i.e. 0.04MJ^Jt, is massive enough to end their lives 
with Type II supernova explosion that peaks at the particle age 
of 5 Myr. The rest of the stellar mass, 0.2 M^q, is considered 
as long-lived stellar mass, M*. [Stage III] At the end of a 
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TABLE 1 
Simulation Suite Description 



Physics" 


MC-TF 


MC-RTF 


Star-forming molecular could (SFMC) particle formation 


(See jl2t 








SFMC feedback: supernova explosion 


(See 323) 








SFMC feedback: ionizing radiation 


(See^Ol 


X 






''For detailed explanation, see the referenced section, o = included, x = not included. 

lifetime of a SFMC particle, 0.76 + 0.04 = 0.8 MJ^jj^ is blown 
back into the ISM with the energy of exploding supernovae. 
Overall, the mass of a SFMC particle is expressed as 



MMc(f)=^Sc 1-0.8 



-&(T'-5)- 



dT' 



1-0. 



0.5-0.5erf{\/8(5-r)} 



(3) 



,(4) 



where erf () is the Gauss error function, and T — t — t^^ is 
the particle age in Myr with particle creation time t^^. This 
formulation replicates the gradual ejection by supernovae of 
80% of MJJjJi that does not end up locked into stars (see Figure 
|2]i. This mass, 0.8 Mj^jt, along with 7.5 x 10"^ ofthe rest mass 
energy of {T ~ 6Myr) — Q.IM^^^ is returned to the cell in 
which the particle resides (see the right panel of Figure [T]i.^ 
The rate of thermal energy input by supernova is proportional 

to which peaks at T 5 Myr with a width of ^ 1 

Myr. Lastly, 2% of the ejected mass is considered as metals, 
contributing to the ISM enrichment. 

To summarize, we emphasize that our model incorporates 
three key observations of star formation in GMCs: (a) the 
slow and inefficient star formation in molecular clouds, (b) 
the commencing of Type II supernova at the death of massive 
stars, and (c) the recycling of a large fraction of the cloud gas, 
blown back into the ISM. 

2.4. Ionizing Radiation Feedback 

Heating by photoionizing radiation from young massive 
stars markedly alters the environment of stellar nursery, and 
it has long been postulated as one of major drive rs of self- 
regulated galactic star formation (IW hitworth ll979l: Williams 
& McKee 1997l:lMatznedl2002h . In this work, we perform an 
explicit, three dimensional transport calculation of ionizing 
radiation from SFMC particles in order to evolve the radiation 
fields throughout the galactic ISM . While the radiative trans- 
fer mac hinery is well described in previous work (e.g. A bel & 
Wandelt' l2002HWise & Abelll2008Ll201 UlKim et al.l20Tli) . we 
briefly describe the basic details relevant to this study. 

Once a SFMC particle is created, ionizing radiation lumi- 
nosity from the particle / is assigned by 



(5) 



where ^mc = 6.3 x lO'*^ photons s 'Mq' is the lifetime- 
averaged ionizing luminosity per solar mass in clusters 

* Our choice of an IMF reaching 300 M ; i instead of the more usual choice 
of 120 M . , is motivated by recent observations suggesting that the IMF is not 
truncated at 100-150 M. (Crowther et al. 2010), as had been claimed based 
on earlier observations. 

' lO"" ergs per every 750 M .i of long-lived stellar mass formed. We adopt 
this rather low value to make supernovae marginally effective in halting the 
runaway collapse of star-forming gas. This choice is designed to further con- 
trast the effect of stellar radiation when two fiducial simulations. MC-TF and 
MC-RTF, are compared. 



(Eq.(lO) of iMurrav & Rahman||2010|) £ph = 16.0eV is flie 
mean energy per dep osited photon (^ Whalen et al.l 120041 : 
IWhalen & Normanl l2006) fixed for monochromatic rays, and 
Mmc(',0 is the SFMC particle mass as a function of particle 
age (see Figure |2]i. Readers should note that we choose to 
employ My\£{i,t) rather than the stellar mass inside it. Our 
choice gives the effective ionizing luminosity per stellar mass 
of 1.9 - 2.6 X 10"^^ photons s^^Mg' at the birth of a SFMC 
particle, depending on the uncertainty of the lower limit of 
the Salpeter IMF (0.1 - 1.0 M©). While the value is slightly 
larger than what was adopted in some previous studies (e.g. 
^ 10 '*^ photons s'^M,:;:' in 'Fuiita et al.1 120031: iDale et all 
I2OI2I) , it helps to approximately encapsulate various other 
forms of fee dback beyond photoion ization, such as protostel- 
lar ou tflows (iNakamura & Lill2008l), stellar winds dOev et all 
I2OOII) . a nd radiation pressu re (iHopkins et al.ll20TTI) . We also 
note that iDale et aTl (1201 2 |) finds that a factor of two uncer- 
tainties in ionizing luminosity are less crucial in determining 
the amount of ionized gas than the gravitational potential in 
star-forming clouds. We will come back to this issue in i j5.4l 

The luminosities Lmc(',0 assigned only to SFMC par- 
ticles that are more than 2.5 kpc away from the galactic center. 
This is to exclude the densest portion of stellar disk so that the 
radiative transport calculation is expedited. This way, we also 
concentrate on the evolution of SFMC particles in the galactic 
spiral arms and outer disk. If the galacto-centric distance to 
the particle is greater than 2.5 kpc, then 12 x 4^ rays (Healpix 
level 3) are isotropically cast from the particle. Each ray is 
traced until most of its photons are absorbed or until it reaches 
the edge of the computational domain, while being adaptively 
split into child rays whenever the angular resolution associ- 
ated with it grows larger than a threshold.^ Here, achieving 
high resolution around SFMC particles is critical in order to 
adequately describe the capture of stellar radiation energy. 

From r = to 6 Myrs after the birth of a SFMC particle, 
the emitted photons interact with the interstellar gas (see the 
middle panel of Figure [TJ. This ionizing radiation feedback 
covers the period between the birth of a SFMC particle and 
the Type II supernova explosion by massive s tars in the clus- 
ter (|HoEkinreFin [2001 [sHflionel^ Though both 
hydrogen and heUum are tracked in the chemistry module, 
photons interact only with hydrogen in the ISM via follow- 
ing three mechanisms. First, a ray loses its photons when it 
ionizes hydrogen at the photoionization rate of 



Pin(l-g'"H) 



(6) 



where fin is the number of incoming photons, Th = nnOyidl 

^ A ray is split as soon as the area resolved by the ray becomes larger than 
0.2 {Ax)- of a local cell. However, to reduce the cost of radiation calculation 
the ray is no longer spht when it is more than 5 kpc away from the source, a 
distance far enough for a dwarf-sized galaxy. In another attempt to speed up 
the calculation two rays may merge to form a single ray. Merging of the rays 
occurs when the distance from the emitting source is more than 10 times the 
separation between the two sources. 
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is the optical depth of a cell, hh is the neutral hydrogen num- 
ber density, C7h is the energy-dependent hydrogen photoion- 
ization cross-section dVerner et alj|1996 ). and dl is the path 
length through the cell. The radiation timestep dt^h set by 
the code is typically comparable to the finest hydrodynami- 
cal timestep. Second, the excess energy above the ionization 
threshold heats the cell with the photoheating rate of 



Pin(l-e-"H)(£ph-£„H) 
nH(Aji:)2t/fph 



(7) 



where h = 13.6eV is the ionization threshold. Third, when 
photons are absorbed by the gas cell, the added acceleration 
by radiation pressure onto the cell is calculated as 



aph 



lost-'-'ph 



mcellt/fph Pgas(Ax)3ct/fph 



(8) 



where t/pph is the photon momentum exerted onto the cell in 
lifph, Piost is the number of photons lost in the cell, and f is 
the directional unit vector of the ray. 

It is worth to briefly note a few points on stellar radiation 
feedback before we proceed: (a) Only hydrogen photoion- 
ization and photoheating are taken into account in the re- 
ported simulation and subsequent analysis, (b) Dust absorp- 
tion is not included in our calculation. This would not signifi- 
cantly change our es cape fraction analysis in f|5]as Figure 1 of 
iGnedin et al.l ( l2008l) demonstrates that the absorption by dust 
does not substantially alter the escape fraction of ionizing ra- 
diation above Lyman limit, fcj By the same token, infrared 
photons re-emitted by UV-irradiated dust grains and the mo- 
mentum imparted by them are not explicitly considered in our 
calculation. Partly because of this reason, radiation pressure 
purely by the hydrogen-ionizing photons does not play a ma- 
jor role in reported results. Including the trapped radiation 
pressure by these absorbed and re-scattered IR photons would 
have furth er enhance d the negative feed back effect (e.g . Hop- 



kins et al. 201 lb. R ecent radiation hydrodynamics simula- 



tions dKrumholz & Thompsoni2012 a.b) however suggest that 
the enhancement in radiation force due to trap p ing is mild, 
significantly less than assumed bv lHopkins et al.1 ( 1201 ll) . 

Table [T] summarizes the feedback models we employ in the 
study. In MC-TF, only supernova feedback is included. In 
MC-RTF, both channels of stellar feedback, ionizing radiation 
and supernova explosion, are simultaneously considered. 

3. INITIAL CONDITIONS 

We simulate the evolution of SFMC particles and their em- 
bedding dwarf galaxy at 3.8 pc resolution. The initial condi- 
tion of our experiment is detailed in this section. 

3.1. Data Conversion Pipeline 

The Gadget-2 (ISpringel|l2005h initial condition generator 
StarScreani^ is used to create an A^-body dataset of an ideal- 
ized isolated galaxy, to which we add gas particles by splitting 
collisionless particles. Then, the Gadget-to-Enzo converter 
hullMethod employs IDL functions qhull (Delaunay tessel- 
lation) and qgridS (linear interpolation) to map the particle 
data onto an adaptively refined, oct-tree structured mesh. For 
details of the data conv ersion pipeline we refer the readers to 
iKimet al.l(l2008ll2009h . 

' The authors are grateful to Jay Bilhngs for making his code public and 
available for this study: http://code.google.eom/p/starscream/. 



3.2. Setting Up A Dwarf Galaxy 

We construct a dark matter halo of 2.3 x lO'^M© that fol- 
lows the Navarro-Fr enk- White profile with a concentration 
parameter c = 10 (N avarro et al.l[T997h . Gas grids are gen- 
erated by splitting the particles with an initial baryonic frac- 
tion of 10%, and with initial metallicity Zjnit = 0.003 Z© ev- 
erywhere. In addition, a collective rotation of spin parameter 
A = 0.05 is provided to both gas and dark matter. We then 
produce a star-forming disk galaxy embedded in a dark mat- 
ter halo by performing a coarse resolution (15.2 pc) simula- 
tion for ^ 1 Gyr in a 32^^ root grid box of 1 Mpc^. While the 
galaxy itself occupies only a small portion of this root level 
box, the large box size is selected so that any boundary effect 
is marginalized. When a relaxed, well-defined galactic disk 
has emerged, ^ 1.3 x IO^^Mq is in 3.4 x 10^ star cluster par- 
ticles formed with a coarse-resolution refinement strategy and 
SFMC particle formation recipe. ^ 6.5 x lO^Mp; is in gaseous 
form either in the ISM or in the halo. These are all embedded 
in a halo of 1 .3 x 10^ dark matter particles. 

We now employ high-resolution (3.8 pc) refinement criteria 
to resolve the galaxy down to the size of an individual molec- 
ular cloud, and make the entire galaxy evolve for another ^ 
30 Myrs. Such high resolution is obtained by applying 13 ad- 
ditional levels of refinement along with the formation criteria 
of - 1000 Mg, SFMC particles (see §2.2). By applying high 
resolution refinement only well into the galactic evolution, 
one can save the computational expense of performing an ex- 
tremely high resolution calculation for a galactic dynamical 
time (> 300 Myrs), but can still observe the behavior of star- 
forming gas clumps for their typical lifetime (~ 30 Myrs). 

4. RESULTS I: OVERVIEW 

A suite of simulations with different stellar feedback is per- 
formed to examine the behavior of SFMC particles and their 
embedding galaxy in one self-consistent framework at high 
spatial resolution. We first overview the performed runs. 

4.1. Summary of Performed Runs 

Each of the calculations is carried out on 16 to 32 proces- 
sors of the Ranger cluster'" at the Texas Advanced Comput- 
ing Center at the University of Texas at Austin. In a typical 
simulation, a galaxy is resolved with ~ 2.0 x 10^ total com- 
putational elements (^ 1.7 x 10^ particles and ^ 160^ cells 
in ~ 10^ grids). In a simulation with both channels of stellar 
feedback (MC-RTF in Table[T]), there are about 3000 to 10000 
active SFMC particles that participate in radiative transfer cal- 
culation. It takes ~ 2000 CPU hours on average to evolve this 
system for 1 Myr A portion of the post-production proce- 
dure and data analysis is conducted on the Orange cluster at 
the Kavli Institute for Particle Astrophysics and Cosmology 
at Stanford University. 

4.2. Self- regulated Star Formation 

In order to exclude the period in which the galaxy is trying 
to reach a new equilibrium with freshly imposed high resolu- 
tion, we investigate the galaxy from 10 to 30 Myr into the evo- 
lution in the subsequent analysis. During this period of ^ 20 
Myrs, the run with both stellar radiation and supernova feed- 
back forms 22.4% less steflar mass (MC-RTF 3.57 /Wgyr"' 

Sun Constellation Linux Cluster, Infiniband-connected AMD, 16 cores 
per node, 2 GB memory per core 



6 



J. Kim et al. 




Fig. 3. — Left columns: density-weighted projection of density in the MC-RTF run in a 20 kpc box at 13.3 Myr into the high-resolution evolution (top: face-on 
/ bottom: edge-on). High-density clumps which bear a number of SFMC particles are noticeable. Right columns: density-weighted projection of photoheating 
rates, in which the shaded region along the right side of the galactic plane is pronounced (black dotted box). Note that the photoheating rates are low in the 
galactic center (< 2.5 kpc in radius, dotted circles) because the particles there are not ray-emitting. 



on average) than the run with only supernova feedback (MC- 
TF, 4.60 Moyr^'). Our realistic description of stellar feed- 
back helps to self -regulate galactic star formation, and provide 
a reasonable star formation rate during the analyzed period of 
the simulation. The suppression ratio could have been larger 
had the SFMC particles in the inner disk also radiated ionizing 
photons (< 2.5 kpc from the galactic center; see i ]2.4l ). Star 
formation relations in connection with t he local Ke nnicutt- 
Schmidt relation dKennicutt'et al.ll2007l : iBigiel et al.i 120081) 
will be extensivel y discussed in an accompanying future study 
(iKim et alJl20T2h . For more discussion on self-regulated star 
formation relations, we refer interested readers to this paper. 

5. RESULTS II: ESCAPE OF IONIZING PHOTONS 

As a first application of our sophisticated feedback model, 
we investigate the escape of ionizing stellar radiation in the 
MC-RTF run in which both ionizing radiation and supernova 
feedback are considered (see TablefTli. We start this section by 
reviewing the recent theoretical advances made in the field. 

5.1. Theoretical Background 

Hydrogen ionizing radiation from dwarf-sized galaxies is 
thought to b e an important driver of reionization in early uni- 
verse (e.g. IVIadau et al. 1999; Faucher-Giguere et al.! 20081 ; 
iBunker et al. 2010.; Kuhlen & Faucher-Giguere .201 2h . Its 
contribution may have increased in time since z ~ 6, relative 



to that of quasars and other sources. Also, galactic ultravio- 
let radiation filtered through the intergalactic medium (IGM) 
is a major source of the metagalactic ultraviolet background. 
This background radiation reduces the gas accretion onto less 
massiv e halos and evaporates their existing gas (e.g. Haardt & 
Madau lT996ll2001blCalverlevetal.ll2011h . 

As such, escape of ionizing radiation from a dwarf galaxy 
has been a fiercely researched topic in extragalactic astron- 
omy. However, the estimates of galactic escape fraction vary 
widely, indicating the difficulty to comprehend the behav- 
ior of ionizing photons within a morphologically complicated 
galactic structure. Various observational studies place the 
galactic escape fraction, /esc, at a few percen t for low-z galax- 
ies (e.g. Heckman et al. 2001; Bergvall et al. 2006; Pellegrini 
et al. ,2012) and high-z galaxies (e.g. Giallongo et al. 200^; 
IShaplev et all 120061; llwata et al.l l2b09; Boutsia et al. 2011), 
with suggested decrease in time likely due to more clumpy 
structures in more evolved galaxies (e.g. Inoue et al. 2006). 

Theoretical investigations have also produced a similar 
value for galactic escape fraction. Analytic estimations em- 
ploy idealized galaxy models and/or Monte Carlo radiative 
transport calculation to evaluate the photon escape fraction 
with varying galactic parameters (e.g. [W ood & Loeb 200(3; 
IClarke & O ev 2002; iFernandez & Shuilll201in . Other au- 
thors have acquired the escape fraction by separately post- 
processing a simulated galaxy which did not integrate the 
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Fig. 4. — Left: all-sky map of projected neutral hydrogen number density (Mollweide equal-area projection) from a radiating SFMC particle in the MC-RTF 
run, seen at 100 kpc at 29.9 Myr into the high-resolution evolution. It displays two notable structures: the galactic disk across the mid-plane, and the gaseous 
clump close to this particle in the bottom right. Middle: all-sky map of density-weighted average of neutral hydrogen fraction seen at 100 kpc. Right: all-sky 
map of photon escape fraction at 100 kpc along different lines of sight, /esc ((' , 6 , <P) . Ionizing photons escape neither from the galactic plane nor from the nearby 
gas clump because of their high column densities of neutral hydrogen. 



effect of ionizi ng stellar radiation in the equation of hydro- 
dvna mics (e.g. iFujita et alj l2003t Raz oumov & Sommer- 
Lars en I2007L 120101: lYaiima et al.ll2009L lIoT Tt Paardekooper 
et al. l2011h . While self-consistent radiation hydrodynamics 
simulations haye been used to stu dy reionization and galaxy 
formation at low resolutio n (e.g. Petkoya & Springelll201 It 
iHasegawa & Semehnl2012h . only occasionally haye they been 
adopted to probe galactic escape fra ction with high spatial res- 
olution (e.g. Gnedin et al. 2008; Wi se & Cen|[2009h . 

As described in f|2] our work incorporates hydrodynamics, 
chemistry, star formation, and ionizing radiation from star- 
forming clouds in one self-consistent framework with an un- 
precedented resolution. We thereby do not neglect the hy- 
drodynamic and chemical influence of the radiation on the 
ISM. Compared to preyious studies, the ray-emitting sources 
in our simulations represent < 10000 radiating SFMC parti- 
cles spatially resolyed in 3.8 pc resolution, scattered on a 
galactic disk. Eyolutionary stages during the lifetime of the 
molecular clouds are also resolyed by the corresponding tem- 
poral resolution (see i)2.3| and i)2.4| i. Additionally, our study 
explores a more massiye and well-defined galactic disk than 
the 3 X 10'' to 3 X lO^M© halos considered in IWise & Ceiil 
(|2009|) . therefore complementing their work. Last but not 
least, our newly-deyeloped approach to post-process the sim- 
ulated data facilitates groundbreaking studies such as exam- 
ining the eyolution of escape fraction per star-forming clump. 
In the remainder of ^j5] we focus on the escape of ionizing 
photons and its angular, temporal, and spatial yariations in 
the MC-RTF run, unless stated otherwise. 

5.2. Escape Fraction via All-Sky Projection 

Figure [3] compiles the projected density and photoheating 
rates in the MC-RTF run at 13.3 Myr into the eyolution with 
high-resolution. It is straightforward to locate the actiyely ra- 
diating SFMC particles around which the photoheating rates 
peak. Also noticeable is the preferential escape of ionizing 
radiation towards polar directions. In particular, the shaded 
region along the right side of the galactic disk plane is highly 
pronounced in the bottom right image of the Figure [3] (black 
dotted box). Note that the photoheating rates are relatiyely 
low in the galactic center because the particles there are not 
emitting radiation (black dotted circles, < 2.5 kpc from the 
galactic center; see ^3.A\ . Radiation fields such as photoheat- 
ing rates (summed for all incoming photons to the cell) are 
tracked by the radiation solyer in each computational cell, and 
later used in chemistry and cooling calculations. 

We, howeyer, aim to eyaluate the fraction of ionizing pho- 
tons that escapes from an indiyidual star-forming clump, or 
from a galaxy. While informatiye, integrated radiation fields 



such as in Figure |3] are less illuminating for our purpose. To 
compute the escape fraction of photons from an indiyidual 
SFMC particle, all-sky maps of the following quantities cen- 
tered on a SFMC particle are generated and shown in Figure 
m (a) the column density of the number of neutral hydro- 
gen, E„jj.iookpc, (b) the density-weighted ayerage of neutral 
hydrogen fraction, and (c) the photon escape fraction along 
different lines of sight. Each of these maps are seen at 100 
kpc away from the ray-emitting source, and at 29.9 Myr into 
the high-resolution eyolution. The angular resolution of these 
maps, A^res = 12x4^ = 12x1 6^ = 3072, corresponds to the 
Healpix leyel — 4 or A^side = 16. The leftmost panel of the pro- 
jected number density shows two prominent gas structures: 
the galactic disk across the mid-plane, and the gaseous clump 
in the bottom right of the map that happens to be yery close 
to this particular source. 

Most notably, the rightmost panel of Figure|4]illustrates the 
all-sky map of escape fraction at 100 kpc along different lines 
of sight from the particle. That is 

/esc(/,0,(p)=e""""""'P^''''^''P' (9) 
where / is the source identification number, d and 9 are the 
angles to the pixel, TH.iookpc = crHJ^»H,iookpc is the optical 
depth at 100 kpc. It is eyident that ionizing photons are pro- 
hibited from escaping the galactic plane or the nearby gas 
clump due to their high column densities. Because the ini- 
tial rays are cast isotropically (see ' ^2A\ , the oyerall escape 
fraction of ionizing photons from this particular source is eas- 
ily obtained by ayeraging /esc(', 9, <p) oyer the entire sky, or 
equiyalently, oyer N-ces equal-area pixels. Hereafter the es- 
cape fraction per SFMC particle, fesc{i), is defined as such; 
in other words, a fraction of photons that reached a 100 kpc 
sphere without being absorbed: 

/esc(0 = ^ T'' r f,scii,e,(p)smeded(p (lo) 

4;r Jo Jo 

«T^I/esc(/,0,(p). (11) 

Further, the escape fraction of a galaxy, /esc, is defined as the 
ayerage of /esc (0 oyer the entire set of ray-emitting sources 
weighted by the initial photon number from each source, 
qMcMucii) inEq.©: 



/esc 



52 ?McA/Mc(0/esc(') 



(12) 



Note that simply by integrating the neutral hydrogen num- 
ber density along different lines of sight we haye acquired 
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Fig. 5. — Left: the locations of SFMC particles sized by their escape fractions, /csc(')' colored by their ages. /esc(')'^ ^■'^ evaluated at 30.7 Myr into the 
high-resolution evolution, and plotted on top of a black-and-white face-on projection of gas density in a 20 kpc box. A single-pixel dot means a SFMC particle 
with zero escape fraction. Right: the weighted sum of the fluxes of escaped photon numbers from SFMC particles along different lines of sight, and seen at 100 
kpc from the ray-emitting sources. The latitude corresponds to the direction vertical to the galactic disk plane, while 0° in the longitudinal coordinate correlates 
with the left side of the galactic plane in the left panel (blue arrow in each panel). The galactic escape fraction at this snapshot is /esc = 0.56%. 



TH,iookpc('i and thus /esc (0- Our post-production pro- 

cedure does not include any additional radiative transfer, and 
is interactive and highly user-adaptive. ' ^ However, while our 
calculation is done in post-processing, we demonstrate in ^5.6\ 
that the result we obtain depends critically on including the 
photoionizing radiation self-consistently during the main cal- 
culation. 

5.3. Escaped Ionizing Photon Flux from A Galaxy and Its 
Angular Variation 

Displayed on the left panel of Figure |5]are the locations of 
SFMC particles whose marker sizes are proportional to the es- 
cape fractions from the particles, /esc(0' 30.7 Myr into the 
high-resolution evolution. The panel depicts the major con- 
tributors of high /esc (0 to the overall galactic escape fraction, 
/esc ~ 0.56%. The age of each particle is represented by the 
color of the circle. Note that the inner disk is not populated by 
radiating SFMC particles because the particles there are not 
ray-emitting (black dotted circle, < 2.5 kpc from the galactic 
center; see ii2.4l l. The galactic escape fraction is heavily dom- 
inated by only a small number of SFMC particles with high 
/esc(0' while the majority of SFMC particles have negligible 
escape fraction as shown by points with a single pixel. We 
will come back to this in ! j5.5l for more discussion. 

For now we concentrate primarily on the galactic escape 
fraction and the total flux of escaped photons. The right panel 
of Figure|5]demonstrates the weighted sum of the fluxes of es- 
caped photon numbers from all SFMC particles along differ- 
ent lines of sight, measured at 100 kpc from the ray-emitting 
sources. This could be interpreted as a escaped photon flux 
from the galaxy, assuming that 100 kpc distance is far enough 
to ignore the separation among ray-emitting SFMC particles. 
Note that 0° in the longitudinal coordinate corresponds to the 
left side of the galactic disk plane in the left panel (blue ar- 
row in each panel), and ±180° to its right side. Therefore, 

' ' Alternatively one may explicitly count the number of photons that have 
been traced for more than 100 kpc, within the radiative transfer simulation. 
This type of calculation is non-interactive and may add additional computa- 
tional cost. 



the prominent shaded region along the galactic plane between 
<p = —180° and — 90° points to the right side of the galactic 
disk (blue dotted box in each panel). As seen in the blue dot- 
ted box, the right side of the galaxy hosts only one distinct 
star-forming clump with many high /esc(0 SFMC particles. 
Photons emanating from this clump is often blocked by thick 
gas spiraling around it. 

Throughout ^ 20 Myrs of high-resolution evolution, we 
find that the flux of escaped photons is in general the small- 
est along the galactic disk plane. Moreover, the escape frac- 
tion varies by up to 4 to 5 orders of magnitude along differ- 
ent lines of sight. For example. Figure |6] shows the escaped 
ionizing photon fluxes at two different epochs. At 17.7 Myr 
ionizing photons rarely escapes along the disk plane of the 
simulated galaxy. At 19.6 Myr photons do escape towards the 
right side of the galaxy, but not through its left side. Zero 
or negligible escape fraction is found due to the large col- 
umn density of neutral hydrogen in a well-defined galactic 
disk, which is seldom observed in prior numerical studies 
of cosmolo2ical 2alaxies (e.2. Gnedin et al. 2008; Wise & 
Cen |2009|) . Nonetheless, the flux of escaped photons is not 
strongly beamed, but rather manifests a large opening angle of 
more than 60° from the galactic pole. We caution the readers 
that the beaming of the escaped photon flux might have been 
enhanced had we allowed SFMC particles in the inner part of 
the galactic disk to radiate. Photon fluxes from the sources in 
the inner disk might have preferentially contributed towards 
the galactic pole thanks to the disk geometry. 

5.4. Temporal Variation of Galactic Escape Fraction 

As discussed in ^15. 21 angular-averaged escape fractions can 
be obtained by integrating the escaped photon flux over the 
entire sky. Figure |7] exhibits the time evolution of galactic 
escape fraction, /esc, from 13.3 Myr to 30.7 Myr of the high- 
resolution evolution. While the mean value of /esc is 1 . 1 %, it 
fluctuates between 0.08% at 17.7 Myr and 5.9% at 24.4 Myr 
The time-averaged /esc of 1.1% is lower than what was pre- 
viously found in simulated low-mass dwarf galaxies at high 
redshift (M < K^-^Mq; IWise & Cenll2009l) . but it is in fine 
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log[Escaped Photons (s 'sr'cm -^)] 
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Fig. 6. — The escaped ionizing photon fluxes measured at 17.7 (left) and 19.6 Myr {right) into the high-resolution evolution. For each snapshot, the galactic 
escape fraction is f^^c = 0.08% and 1.1%, respectively. Note that at 17.7 Myr ionizing photons rarely escapes along the disk plane of the simulated galaxy. See 
the caption of Figure|5]for a detailed description of the coordinate system. 

with the va lues observed in larger galaxies (M > lO'^M©; 
iGnedin et a l. 2008) or in idealized galaxie s with rotationally 
supported disks tPaardekooper et all 1201 Ih . Also found is a 
mild positive correlation between /esc and the instantaneous 
star formation rate, dM^,/dt (with dt = 0.37 Myr), plotted 
with a dashed line. The large variation of /esc up to more than 
an order of magnitude in just a few Myrs might be related to 
the change in instantaneous star formation rate. These obser- 
vations are consistent with the previous escape fraction stud- 
ies using radiation hvdrodvna mics simulations (e.g. Gnedin 
et al. . 2008;. Wise & Cen.200a) . 

We, however, caution that /esc might have changed if the 
SFMC particles in the inner disk had radiated ionizing pho- 
tons. Other sources of uncertainty include: (a) the change 
in supernova energy that disperses the surroundi ng gas, (b) 
the ch ange in ionizing luminosity, Eq.©, although lDale et all 
argues that a factor of a few uncertainties in ionizing 
luminosity are not of significant importance in determining 
the amount of ionized gas in star-forming clouds, and (c) the 
inclusion of dust absorption, although Gnedin et al.l (l2008h 
claims that it does not substantially change /esc above Lyman 
limit (see ^2A[ . 
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5.5. Escape Fraction Per SFMC Particle 

Two critical factors that distinguish our work from previous 
escape fraction studies are that: faj the ray-emitting sources in 
our simulations model star-forming molecular clouds, and (b) 
our new post-production procedure enables us to measure the 
evolution of escape fraction per SFMC particle. Understand- 
ing how the ionizing stellar radiation feedback manifests itself 
on an evolving star-forming clump in a galactic context is of 
vital interest, especially in comparison with GMC-scale stud- 
ies in which ionizing feedback is considered (e.g. Krumholz 
et al .l2006, ,2007.: .Mac L ow et al. 2007 ;,Vazquez-S emadeni 
et al. 120 lot [Arthur et al.ll20 1 1 : Goldbau metaLllMTl: Dale & 
Bonnell l201 iHDale et al.ll2012i:iGendelev & Krumholzl2012h . 
Moreover, observationally, stellar radiation leaked from in- 
dividual molecular clouds is thought to power the ionization 
of diffuse ionized gas (DIG) in a galaxy, and maintain s DIG 
layer above and below the galactic disk (e.g. Revnoldd l 19841 : 
iHoopes & Walterbosll2003l:IVoges & Walterbosii2006) . In this 
section, we inspect the escape fraction per SFMC particle, 
/esc(0' and how it e volv es as the particle ages. 

As examined in ^ ^5.31 the galactic escape fraction is dom- 
inated by a small number of SFMC particles with very high 
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Time into High-resolution Evolution (iviyr) 

Fig. 7. — Solid line: time evolution of galactic escape fraction, /esc, from 
13.3 Myr to 30.7 Myr of die high-resolution evolution in the MC-RTF run. It 
fluctuates between 0.08% at 17.7 Myr and 5.9% at 24.4 Myr. Dashed line: 
instantaneous star formation rate in the unit of M0yr^ ' . 

fesc{i)- At the same time a large fraction of SFMC particles 
have zero or negligible escape fraction as shown in Figure |5] 
by points with a single pixel. This observation is dramatically 
displayed in Figure [8] where actively radiating SFMC parti- 
cles are scattered on the plane of particle age and the escape 
fraction per particle. This snapshot is made at 29.9 Myr into 
the high-resolution evolution. It is worth noting that out of 
8709 actively radiating SFMC particles at this point only 140 
have /.sc (0 greater than 1 %. This is merely 1 .6% of the entire 
population of SFMC particles. Among these the variation of 
/esc(0 is broad, from just ^ 1% up to 80%, a feature sim- 
ilar to what Dal e et al.l (120121) found. On the other hand, the 
majority of the SFMC particles fall on top of the jc-axis on the 
left panel of Figure[8]with negligible /esc(0- This observation 
strongly supports the argument that the galactic escape frac- 
tion is determined by a few SFMC particles with exception- 
ally high escape fractions. (The following toy model may help 
elucidate our finding: if galactic /esc is 1%, then it is more 
plausible that 1% of the SFMC particles have /esc(0 = 100%, 
rather than 100% have /esc(0 = !%■) From this observation, 
one may claim that the escape of ionizing photons is primarily 
determined by the small-scale properties of the star-forming 
clouds in the vicinity of a SFMC particle, not by the proper- 
ties of the galaxy on larger scales. One may also argue that 
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Fig. 8. — The distribution of actively radiating SFMC particles (i.e. particle age T(i) < 6 Myrs) on the plane of the particle age and the escape fraction 
per particle, /esc(0> nieasured at 29.9 Myr into the high-resolution evolution with two types of y-axis, linear (left) and logarithmic (right). At this snapshot the 
galactic escape fraction is /esc = 0.61%. Note that most of the SFMC particles fall on the .v-axis on the left panel with negligible fesci')', only 1.6% of the actively 
radiating particles are shown to have /csc(0 greater than 1%. Note also that the escape fraction slowly rises as a function of particle age. This trend is in part 
aided by the supernova explosion feedback that peaks at T{i) = 5 Myr that is marked with a dotted vertical line. 



galactic escape fraction cannot be properly calculated without 
self-consistently considering the interaction of the ionizing ra- 
diation at molecular cloud scales. 

Sometimes several SFMC particles cluster in or around a 
single computational cell and show nearly identical escape 
fractions. This phenomenon is manifested in Figure |8] where 
18 particles have /esc(0 ^ 80%. This type of cell has ended 
up hosting multiple SFMC particles because it sits right at the 
center of a gravitationally unstable gas clump that potentially 
exhibits runaway collapse. This cluster of SFMC particles can 
be regarded as a large stellar nursery or a massive young star 
cluster of > lO^M© that is under-resolved with our numerical 
resolution, 3.8 pc. 

Readers should also note that the escape fraction slowly 
rises as a function of particle age. By averaging 20 snapshots 
between 13.3 Myr and 30.7 Myr into the high-resolution evo- 
lution, we discover that /esc(0 increases from 0.27% at its 
birth to 2.1% at the end of a SFMC particle lifetime, 6 Myrs. 
In other words, our best fit is 

/esc(0 = 0.0031 X r(/) +0.0027 (13) 

where T{i) — t ~ fcr(/) is the age of the particle in Myr with 
particle creation time fcr(0' defined in ^ ^2.31 This result in- 
dicates that young SFMC particles are usually still buried in 
cold neutral gas in which they were spawned. As the SFMC 
particles age, they may drift away from the dense gas clump 
causing /esc (0 to rise. The trend is also aided by ionizing stel- 
lar radiation and supernova explosion feedback (peaking at 5 
Myr marked with a vertical line in Figure |8]l which disperse 
the gas around the particle. 

Because older SFMC particles tend to have higher /esc(0' 
they contribute equally or more to the total escaped photons 
despite the fact that they radiate fewer photons at the outset. 
Note that Mmc(' ,0 in Eq.© is a function of time decreasing 
sharply from Mjjjl at T = 4 Myr to 0.2 Mjjjt at T = 6 Myr. 
The radiation luminosity, Lmc , ? )' experiences the same 80% 
reduction in the last 2 Myrs of SFMC particle lifetime. How- 
ever the increase in /esc(0 typically offsets the decrease in 
Lmc('iO' so the number of escaped photons scarcely depend 
on the particle age. In order to illustrate this discovery, cu- 
mulative numbers of escaped ionizing photons are drawn as a 



function of /esc (0 in Figure|9] For all SFMC particles at 29.9 
Myr into the high-resolution evolution, we compute both the 
number of ionizing photons which escape from that particle, 
and /esc(0- We then break the particles into three age bins of 
[0, 2] Myr, [2, 4] Myr, and [4, 6] Myr. Within each bin, we 
order the particles by escape fraction, from lowest to highest, 
then compute how many escaping photons come from parti- 
cles with that value of /esc (0 or less, and plot the result in 
Figure |9] At this particular timestep, the contribution to the 
total escaped photons by the particles of age [4, 6] Myr is the 
greatest, close to a half of the total number It is followed by 
the other age bins, [0, 2] Myr and [2, 4] Myr Again, this 
phenomenon is in spite of the fact that the particles of age [4, 
6] Myr originally radiate the least number of photons. Our 
result bears a striking rese mblance to t he rad iation hydrody- 
namics simulation result by lDale et al.l (1201 2h . where clusters 
with the lowest initial luminosities have the highest escape 
fractions, and vice versa. These observations suggest that the 
escaped stellar photons may at least provide a sufficient num- 
ber of photons to sustain galactic DIG layers. 

Lastly, various other useful information can also be ex- 
tracted by carefully examining this figure. For example, for 
the age bin of [0, 2] Myr, particles with /esc (0 ^ 80% con- 
tributes the most to the escaped photons in this bin. In con- 
trast, the contribution by particles with /esc(') < 75% is less 
than 1%. In all three age bins, contributions by particles with 
/esc(0 < 20% is no more than ^ 5%. This assessment reaf- 
firms the claim that the escaped photon flux from a galaxy is 
predominantly controlled by a small number of SFMC parti- 
cles with exceptionally high /escfO- It also is consi stent with 
the observational findings bv lPellegrini et al.l (l2012l Fig. 21). 

5.6. Comparison to the Run Without Stellar Radiation 

We have thus far focused on the MC-RTF run, which in- 
cludes stellar radiation feedback. However, it is very illumi- 
nating to compare the results to the MC-TF run, which in- 
cludes supernova feedback but not radiation. To examine this, 
we have post-processed this run in exactly the same manner 
as the MC-RTF run. Surprisingly, the result of this exercise is 
that the escape fraction of the MC-TF run is negligible at all 
timesteps. Figure[TO]shows an example, comparing runs MC- 
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0.4 0.6 
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Fig. 9. — The cumulative numbers of escaped ionizing photons as a function of escape fraction per SFMC particle, fssd')- Derived at 29.9 Myr into the 
high-resolution evolution. Thi'ee different age bins are shown. For example, for the age bin of [0, 2] Myr, pailicles with /esc Q) ^ 80% contributes the most to the 
escaped photons in this bin. In all three age bins, contributions by particles with /c.sc(') < 20% are minimal. More examination of this plot is provided in the text. 



TF and MC-RTF at 30.7 Myr into the high-resolution evolu- 
tion. This difference is not because the behaviors of star for- 
mation in the two runs are radically different. Instead, it is be- 
cause supernova feedback alone fails to provide sufficient ion- 
ization to create paths of low neutral fraction through which 
photons can escape. This is partly a matter of timing: by the 
time supernovae going off and clearing the gas around SFMC 
particles, the ionizing luminosity is already in decline. It is 
also partly a matter of ionization budget: supernovae carry 
a lot of energy, but they tend to produce small amounts of 
hot gas but leave large quantities of neutral gas around them, 
and the neutral gas effectively blocks the escape of ionizing 
photons. The comparison of runs MC-RTF and MC-TF illus- 
trates that one can only capture the escape fraction properly 
by self-consistently including photoionization as part of the 
calculation. Post-processing alone is insufficient. 

6. DISCUSSION 

6.1. Summary and Conclusions 

We have performed a comprehensive high-resolution sim- 
ulation of a dwarf galaxy including a sophisticated model of 
stellar feedback. Our goal has been to investigate galactic star 
formation and the formation of a dwarf-sized galaxy in a sin- 
gle self-consistent framework, without omitting the detailed 
interaction of star-forming gas clumps with their surrounding 
galactic envelope. Our major findings are as follows. 

1 . Self- consistent Numerical Study of Galactic Star For- 
mation: We have implemented a new way of describing 
feedback from SFMC particles in galactic simulations 
by combining ionizing stellar radiation and supernova 
explosion. The ionizing radiation feedback from each 
of < 10000 SFMC particles is rendered by tracing the 
ultraviolet photon rays on the fly (f|2]l. Joined with high 
numerical resolution of 3.8 pc, the realistic description 
o f ste llar feedback helps to self-regulate star formation 
( ii4.2b . Our simulation incorporates galactic gas enve- 
lope and SFMC particles in one numerical framework, 
and permits us to explore how galaxies and star clusters 
co-evolve and influence each other. 

2. Escape Fraction via All-Sky Projection: By integrat- 
ing the neutral hydrogen number density along different 



lines of sight we obtain the escape fraction of ionizing 
radiation from an individual SFMC particle or from a 
galaxy ( i)5.2t . Our post-production procedure does not 
include any additional radiative transfer, and is interac- 
tive and highly user-adaptive. The escape fraction study 
in this work is possible only because three key compo- 
nents are consistently integrated with one another: (a) 
the realistic stellar feedback of ionizing radiation, (b) 
high spatial resolution and primordial chemistry to re- 
produce the galactic ISM, and (c) the post-production 
procedure to estimate the flux of escaped photons from 
an individual SFMC particle, or from a galaxy. The ma- 
chinery established in this work brings a unique insight 
into the physics of photon escape fraction. 

3. Angular and Temporal Variations of Galactic Escape 
Fraction: By simulating a galactic halo of 2.3 x 

we find that the galactic escape fraction, /esc, 
fluctuates between 0.08% and 5.9% during a 20 Myr 
period with a mean value of 1 . 1 % ( ij5.4t . The flux of es- 
caped photons is not strongly beamed, but manifests a 
large opening angle of more than 60° from the galactic 
pole. It is often the smallest along the galactic dis k du e 
to the large column density of neutral hydrogen (^53]l. 

4. Escape Fraction Per Star-forming Particle: We also in- 
spect the escape fraction per SFMC particle, /esc (')■ We 
discover that the galactic escape fraction is dominated 
by a small number of SFMC particles with exception- 
ally high /esc(0- Oils i^^y therefore argue that the es- 
cape of ionizing photons is primarily determined by the 
small-scale properties of the clouds in the vicinity of a 
young star cluster, not by an overall galactic structure. 
On average, the escape fraction from a SFMC particle 
rises from 0.27% at its birth to 2.1% at the end of a 
particle lifetime, 6 Myrs. This is because SFMC par- 
ticles drift away from the dense gas clumps in which 
they were born, and because the gas around the SFMC 
particles is dispersed by ionizing radiation and super- 
nova feedback. Since older SFMC particles tend to 
have higher fesc{i), they contribute equally or more to 
the total escaped photons despite the fact that they ini- 
tially radiate fewer photons ( ^5.5l l. 
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Fig. 10. — The escaped ionizing photon fluxes measured at 30.7 Myr into the high-resolution evolution for the MC-TF (left, only with supernova feedback) and 
MC-RTF runs (right, with both supernova and stellar radiation feedback). The galactic escape fraction is /^sc = 0.00% and 0.56%, respectively. See the caption 
of Figure|5]for a detailed information about this projection. For the MC-TF run the map shows what the photon fluxes would have looked like had the SFMC 
particles begun to radiate at this very moment; note that for this run, fluxes in all pixels are considerably lower than the lower bound of the color bar. Because 
SFMC particles are often buried deep in the cold neutral gas, photons would have rarely escaped the galaxy. 



6.2. Future Work 

There is a large, multi-dimensional parameter space to ex- 
plore in connection with the galactic and clump escape frac- 
tion: spatial resolution, mass of a galactic halo, strength of 
stellar feedback , drift of runaway young stars (e.g. Ceverino 
et al.^010), just to name a few. The new feedback scheme 
described in this study also opens numerous doors not only to 
understand the escape of ionizing photons from star-forming 
clumps, but to make mock obseryation such as Ha emis- 
sion and to examine eyolying enyironment of star-forming 
gas clumps. Based on the physics and framework deyeloped 
in this work, galactic star formation and star formatio n rela- 
tions are studied in the future study dKim et al.ll2012l) . The 
accompanying future paper discusses an analysis of star for- 
mation relations in dwarf-sized galaxies, and how well they 
agree with obseryations. 



J. K. thanks Jeff Oishi, Joel Primack, and Chao-chin Yang 
for proyiding insightful comments and yaluable adyice. M. 



Abel, T, & Wandelt, B. D. 2002, MNRAS, 330, L53 

Anninos, P., Zhang, Y., Abel, T, & Norman, M. L. 1997, New Astronomy, 2, 
209 

Anninos, W. Y., & Norman, M. J. 1994, ApJ, 429, 434 

Arthur, S. J., Henney, W. J., Mellema, G., de Colle, F.. & Vazquez-Semadeni. 

E. 2011, MNRAS, 414, 1747 
Bergvall, N., Zackrisson, E., Andersson, B.-G., Amberg, D., Masegosa, J., & 

Ostlin, G. 2006, A&A, 448, 513 
Bigiel, F, Leroy, A., Walter, F, Brinks, E., de Blok, W. J. G., Madore, B., & 

Thornley, M. D. 2008, AJ, 136, 2846 
Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F, & Blitz, L. 2008, 

ApJ, 686, 948 

Boumaud, F, Elmegreen, B. G., Teyssier, R., Block, D. L., & Puerari, 1. 2010, 
MNRAS, 409, 1088 

Boutsia, K., Grazian, A., Giallongo, E., Fontana, A., Pentericci, L., 
Castellano, M., Zamorani, G., Mignoh, M., Vanzella, E., Fiore, F, Lilly, 
S. J., Gallozzi, S., Testa, V., Paris, D., & Santini, P 2011, ApJ, 736, 41 

Bryan, G. L., & Norman, M. L. 1997, in Astronomical Society of the Pacific 
Conference Series, Vol. 123, Computational Astrophysics; 12th Kingston 
Meeting on Theoretical Astrophysics, ed. D. A. Clarke & M. J. West, 363- 
368 

Bunker, A. J., Wilkins, S., Elhs, R. S., Stark, D. P, Lorenzoni, S., Chiu, K., 
Lacy, M., Jarvis, M. J., & Hickey, S. 2010, MNRAS, 409, 855 

Calverley, A. P, Becker, G. D., Haehnelt, M. G., & Bolton, J. S. 2011, 
MNRAS, 412, 2543 

Cen, R., & Ostriker, J. R 1992, ApJ, 399, LI 13 

Ceverino, D., Dekel, A., & Boumaud, F 2010, MNRAS, 404, 2151 

Clarke, C, & Oey, M. S. 2002, MNRAS, 337, 1299 



R. K. acknowledges support from an Alfred P. Sloan Fellow- 
ship, from the NSF through grant CAREER-0955300, and 
from NASA through a Chandra Space Telescope Grant and 
through Astrophysics Theory and Fundamental Physics Grant 
NNX09AK31G. J. H. W. gratefully acknowledges support 
from the NSF Grant AST-1211626. M. J. T. gratefully ac- 
knowledges support from the NSF Grant OCl- 1048505. N. J. 
G. is supported by a Graduate Research Fellowship from the 
NSF. The examination of the simulation data and the post- 
production analysis are grea tly aided by an AMR analysis 
toolkit yt (iTurk et al.ll2011h . This work used the Extreme 
Science and Engineering Discoyery Enyironment (XSEDE), 
which is supported by NSF grant OCl-1053575. The au- 
thors acknowledge the Texas Adyanced Computing Center 
(TACC) at the Uniyersity of Texas at Austin for proyiding 
high-performance computing resources that haye contributed 
to the research results reported within this paper. The authors 
are also grateful for the support from Stuart Marshall, Ken 
Zhou, and the computational team at SLAC National Accel- 
erator Laboratory during their usage of the Orange cluster 



Crowther, P. A., Schnurr, O., Hirschi, R., Yusof, N., Parker, R. J., Goodwin, 

S. R, & Kassim, H. A. 2010, MNRAS, 408, 731 
Dale, J. E., & Bonnell, 1. 2011, MNRAS, 414, 321 
Dale, J. E., Ercolano, B., & Bonnell, 1. A. 2012, MNRAS, 424, 377 
Dobbs, C. L., Burkert, A., & Pringle, J. E. 2011, MNRAS, 413, 2935 
Faucher-Giguere, C.-A., Lidz, A., Hemquist, L., & Zaldarriaga, M. 2008, 

ApJ, 688, 85 

Fernandez, E. R., & Shull, J. M. 201 1, ApJ, 731, 20 

Fujita, A., Martin, C. L., Mac Low, M.-M., & Abel, T. 2003, ApJ, 599, 50 

Fukui, Y, Kawamura, A., Wong, T., Mural, M., Iritani, H., Mizuno, N., 

Mizuno, Y, Onishi, T, Hughes, A., Ott, J., Muller, E., Staveley-Smith, 

L., & Kim, S. 2009, ApJ, 705, 144 
Gendelev, L., & Ki-umholz, M. R. 2012, ApJ, 745, 158 
Giallongo, E., Cristiani, S., D'Odorico, S., & Fontana, A. 2002, ApJ, 568, L9 
Gnedin, N. Y, Kravtsov, A. V., & Chen, H.-W. 2008, ApJ, 672, 765 
Goldbaum, N. J., Ki-umholz, M. R., Matzner, C. D., & McKee, C. F 2011, 

ApJ, 738, 101 

Governato, F, Zolotov, A., Pontzen, A., Christensen, C, Oh, S. H., Brooks, 

A. M., Quinn, T., Shen, S., & Wadsley, J. 2012, MNRAS, 422, 1231 
Haardt, F, & Madau, P 1996, ApJ, 461, 20 

Haardt, F, & Madau, P. 2001, in Clusters of Galaxies and the High Redshift 

Universe Observed in X-rays, ed. D. M. Neumann & J. T. V. Tran 
Hasegawa, K., & Semehn, B. 2012, ArXiv e-prints, 1209, arxiv: 1209.4143 
Heckman, T. M., Sembach, K. R., Meurer, G. R., Leitherer, C, Calzetti, D., 

& Martin, C. L. 2001, ApJ, 558, 56 
Hoopes, C. G., & Walterbos, R. A. M. 2003, ApJ, 586, 902 
Hopkins, P F, Quataert, E., & Murray, N. 201 1, MNRAS, 417, 950 



Dwarf Galaxies with Radiation Feedback. I 



13 



Hughes, A., Wong, T., Ott, J., MuUer, E., Pineda, J. L., Mizuno, Y., 
Bernard, J.-R, Paradis, D., Maddison, S., Reach, W. T., Staveley-Smith, 
L., Kawamura, A., Meixner, M., Kim, S., Onishi, T., Mizuno, N., & Fukui, 
Y. 2010, MNRAS, 406, 2065 
Inoue, A. K., Iwata, I., & Deharveng, J.-M. 2006, MNRAS, 371, LI 
Iwata, I., Inoue, A. K., Matsuda, Y, Furusawa, H., Hayashino, T., Kousai, K., 
Akiyama, M., Yamada, T., Burgarella, D., & Deharveng, J.-M. 2009, ApJ, 
692, 1287 

Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 

Kennicutt, Jr., R. C, Calzetti, D., Walter, P., Helou, G., Hollenbach, D. J., 
Armus, L., Bendo, G., Dale, D. A., Draine, B. T., Engelbracht, C. W., 
Gordon, K. D., Prescott, M. K. M., Regan, M. W., Thornley, M. D., Bot, C., 
Brinks, E., de Blok, E., de Mello, D., Meyer, M., Moustakas, J., Murphy, 
E. J., Sheth, K., & Smith, J. D. T. 2007, ApJ, 671, 333 

Kim, J.-H., Krumholz, M. R., Wise, J. H., Turk, M. J., Goldbaum, N. J., & 
Abel, T. 2012, ApJ, to be submitted 

Kim, J.-H., Wise, J. H., & Abel, T. 2008, in American Institute of Physics 
Conference Series, Vol. 990, First Stars 111, 429^31 

Kim, J.-H., Wise, J. H., & Abel, T. 2009, ApJ, 694, L123 

Kim, J.-H., Wise, J. H., Alvarez, M. A., & Abel, T. 2011, ApJ, 738, 54 

Koyama, H., & Inutsuka, S. 2002, ApJ, 564, L97 

Krumholz, M. R., Matzner, C. D., & McKee, C. F. 2006, ApJ, 653, 361 

Krumholz, M. R., Stone, J. M., & Gardiner, T. A. 2007, ApJ, 671, 518 

Krumholz, M. R., & Tan, J. C. 2007, ApJ, 654, 304 

Krumholz, M. R., & Thompson, T. A. 2012a, ArXiv e-prints, 1203, 

arXiv:1203.2926 
— . 2012b, in preparation 

Kuhlen, M., & Faucher-Gigufere, C.-A. 2012, MNRAS, 423, 862 

Li, Y, Mac Low, M.-M., & Klessen, R. S. 2005, ApJ, 626, 823 

Lombardi, M., Alves, J., & Lada, C. J. 2010, A&A, 519, L7 

Mac Low, M.-M., Toraskar, J., Oishi, J. S., & Abel, T. 2007, ApJ, 668, 980 

Madau, R, Haardt, R, & Rees, M. J. 1999, ApJ, 514, 648 

Matzner, C. D. 2002, ApJ, 566, 302 

Murray, N., Quataert, E., & Thompson, T. A. 2010, ApJ, 709, 191 

Murray, N., & Rahman, M. 2010, ApJ, 709, 424 

Nakamura, R, & Li, Z.-Y. 2008, ApJ, 687, 354 

Navarro, J. R, Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Norman, M. L., & IJryan, G. L. 1999, in Astrophysics and Space Science 
Library, Vol. 240, Numerical Astrophysics, 19 

Norman, M. L., Bryan, G. L., Harkness, R., Bordner, J., Reynolds, D., 
O'Shea, B., & Wagner, R. 2007, ArXiv e-prints, 705, arXiv:0705.1556 

Oey, M. S., Clarke, C. J., & Massey, P. 2001, in Dwarf galaxies and then- 
environment, ed. K. S. de Boer, R.-J. Dettmar, & U. Klein, 181 

Paardekooper, J.-R, Pelupessy, F. I., Altay, G., & Kruip, C. J. H. 2011, A&A, 
530, A87 



Pellegrini, E. W., Oey, M. S., Winkler, R R, Points, S. D., Smith, R. C, 

Jaskot, A. E., & Zastrow, J. 2012, ApJ, 755, 40 
Petkova, M., & Spnngel, V. 201 1, MNRAS, 412, 935 

Plume, R., Jaffe, D. T., Evans, 11, N. J., Martin-Pintado, J., & Gomez- 
Gonzalez, J. 1997, ApJ, 476, 730 
Razoumov, A. O., & Sommer-Larsen, J. 2007, ApJ, 668, 674 
—.2010, ApJ, 710, 1239 
Reynolds, R. J. 1984, ApJ, 282, 191 
Salpeter, E. E. 1955, ApJ, 121, 161 
Schmidt, M. 1959, ApJ, 129, 243 

Shapley, A. E., Steidel, C. C, Pettini, M., Adelberger, K. L., & Erb, D. K. 

2006, ApJ, 651,688 
Shetty, R., & Ostriker, E. C. 2008, ApJ, 684, 978 
Springel, V. 2005, MNRAS, 364, 1105 

Springel, v., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 776 
Stinson, G., Brook, C, Maccio, A. V., Wadsley, J., Quinn, T. R., & 

Couchman, H. M. R 2012, ArXiv e-prints, 1208, arXiv: 1208.0002 
Stone, J. M., & Norman, M. L. 1992a, ApJS, 80, 753 
— . 1992b, ApJS, 80, 791 

Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 

Tasker, E. J., & Tan, J. C. 2009, ApJ, 700, 358 

Teyssier, R., Chapon, D., & Boumaud, F. 2010, ApJ, 720, L149 

Truelove, J. K., Klein, R. I., McKee, C. R, HoUiman, II, J. H., Howell, L. H., 

& Greenough, J. A. 1997, ApJ, 489, L179+ 
Turk, M. J., Smith, B. D., Oi.shi, J. S., Skory, S., Skillman, S. W., Abel, T., & 

Norman, M. L. 201 1, ApJS, 192, 9 
Vazquez-Semadeni, E., Colin, P., Gomez, G. C, Ballesteros-Paredes, J., & 

Watson, A. W. 2010, ApJ, 715, 1302 
Vemer, D. A., Periand, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 

465, 487 

Voges, E. S., & Walterbos, R. A. M. 2006, ApJ, 644, L29 

Whalen, D., Abel, T., & Norman, M. L. 2004, ApJ, 610, 14 

Whalen, D., & Norman, M. L. 2006, ApJS, 162, 281 

Whitworth, A. 1979, MNRAS, 186, 59 

WilUams, J. R, & McKee, C. R 1997, ApJ, 476, 166 

Wise, J. H., & Abel, T. 2008, ApJ, 685, 40 

— . 2011, MNRAS, 414, 3458 

Wise, J. H., & Cen, R. 2009, ApJ, 693, 984 

Wood, K., & Loeb, A. 2000, ApJ, 545, 86 

Yajima, H., Choi, J.-H., & Nagamine, K. 2011, MNRAS, 412, 411 
Yajima, H., Umemura, M., Mori, M., & Nakamoto, T. 2009, MNRAS, 398, 
715 



